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We study the quasi-simultaneous Spectral Energy Distributions (SED) of 48 bright blazars, col- 
lected within three months of the Fermi LAT Bright AGN Sample (LBAS) data taking period, 
combining Fermi and Swift data with radio NIR-Optical and hard-X/7-ray data. We have used 
these SEDs to characterize the peak position and intensity of both the low and the high-energy 
features of blazar spectra. The results have been used to derive empirical relationships that esti- 
mate the position of the two peaks from the broad-band colors (i.e. the radio to optical, Oro, and 
optical to X-ray, Oox, spectral slopes) and from the 7-ray spectral index. Our data show that the 
synchrotron peak frequency is positioned between lO^'^'^ and 10^*'^ Hz in broad-lined FSRQs and 
between 10 13 and 10^^ Hz in featureless BL Lacertae objects. We find that the 7-ray spectral 
slope is strongly correlated with the synchrotron peak energy and with the X-ray spectral index, as 
expected at first order in synchrotron - inverse Compton scenarios. However, simple homogeneous, 
one-zone. Synchrotron Self Compton (SSC) models cannot explain most of our SEDs, especially in 
the case of FSRQs and low energy peaked (LBL) BL Lacs. More complex models involving External 
Compton Radiation or multiple SSC components are required to reproduce the overall SEDs and 
the observed spectral variability. 



INTRODUCTION 



The Fermi Gamma-Ray Space Telescope was 
launched on June 11, 2008 from Cape Kennedy, 
Florida opening a new era in the knowledge of high- 
energy 7-ray band. The high sensitivity, large dinamic 
range and excellent time coverage of the Large Area 
Telescope, the primary instrument on board Fermi, 
all represent significant advantages over the previous 
7-ray observations and will play a key role in the un- 
derstanding of the physical processes underlying the 
high energy emission of blazars. Typical AGN are 
radio-quiet and are found in large numbers of optical 
and X-ray frequencies. Blazars are rare extra-galactic 
objects, about 10% of the entire population of AGN, 
where a jet of plasma, moving at relativistic speed 
is pointing close to the line of sight of the observer 
[l|. However, blazars strong non-thermal emission 
over the entire electro-magnetic spectrum makes them 
the dominant type of extragalactic sources in those en- 
ergy windows where the accretion onto a supermassive 
black hole, or other thermal machenism, do not pro- 
duce much radiation. Although blazars have been ob- 
served for tens of years now, the exististing experimen- 
tal data do not allow one to unambiguously identify 
the physical mechanism responsible for the electro- 
magnetic emission. Given the existing high sensitivity 
detectors to study the low energy component (radio, 



IR, optical and X-rays), one of the main reasons for 
this lack of understanding was the low instrumental 
background between 0.5 MeV to 300 GeV. This poor 
sensitivity often precluded detailed cross-correlation 
studies between low and high energy emission mod- 
els, implying more and higher quality data at 7-ray 
energies to better understand these objects. Simul- 
taneous and detailed Spectral Energy Distributions 
(SEDs) are indeed needed to cover the greatest possi- 
ble band of the electromagnetic spectrum. In fact, one 
of the most effective ways of studying blazar 's physical 
properties is through the use of multi-frequency data 
that, so far, have been carried out almost exclusively 
on the occasion of large flaring events of a few bright 
and well known blazars. 



In this contribution I report on the properties of 
broad-band (radio to high energy 7-ray) SEDs of a 
sample of 48 Fermi bright blazars reported in Q]- 
The detailed SEDs have been built and studied using 
simultaneous or quasi-simultaneous data taken from 
May 2008 to January 2009 and obtained from Fermi, 
Swift, and other space based or ground telescopes such 
as AGILE 7-ray data, Effelsberg, OVRO, Ratan-600 
1-22 GHz, radio, mm, NIR and optical data from the 
GASP-WEBT collaboration, Mid-infrared VISIR ob- 
servations. 
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II. BLAZAR CLASSIFICATION SCHEME 
AND DATA SAMPLE 

The results of Fermi LAT observations in the first 
three months of science operations are reported in the 
first LAT Bright Gamma-Ray Source hst [3] that in- 
cludes 132 sources at galactic latitude |b| >10°, with 
test statistic greater than 100 (> 10 a). Among these, 
a sample of 106 sources indicate high-confidence asso- 
ciations with known AGN and is referred to as the 
LAT Bright AGN Sample (LB AS) 0], comprising 58 
FSRQs, 42 BL Lac objects, 4 blazar of uncertain type 
and 2 radio-galaxies. It is confirmed that the bright 
7-ray sky is dominated by radio-loud AGN. In fact 
about 90% of the LBAS sources is associated to AGN 
listed in radio catalogs (CRATES/CGRaBS, BZCat). 
Moreover only about one third of the bright Fermi 
AGN were also detected by EGRET. 

Blazars represent a subgroup of non-thermal domi- 
nated AGN. For a more detailed classification of AGN 
refer to Q. Blazars are extremely strong and highly 
variable, core dominated flat or inverted spectrum ra- 
dio loud AGN whose jet is looked closely to the line of 
sight of the observer. Blazars can be further divided 
into two subgroups: 

• FSRQs showing broad emission lines in their op- 
tical spectrum just like normal QSO (Quasi Stel- 
lar Objects) and including Flat Spectrum Radio 
Quasars (FSRQs) and broad line radio galaxies. 

• BL Lacs (or BL Lacertae objects) that do not 
show any strong or broad line in their opti- 
cal spectrum and whose radio compactness and 
broad-band SED are very similar to that of 
strong lined blazars. 

A further classification of BL Lac objects com- 
prises three subclasses (LBL IBL and HBL, as Low, 
Intermediate and High peaked BL Lacs) depending 
on the position in the SED of the peak energy of 
the synchrotron emission that reflects the maximum 
energy the particles can be accelerated in the jet. 
This classiflcation has been extended to all blazars 
that can be classifled into Low, Intermediate and 
High energy (Synchrotron) Peaked blazars, respec- 
tively called LSP (synchrotron peak Vp^ak <10^"' Hz), 
ISP (1014 < upeak <10i5) and HSP {vp^ak >IQ^'') In 
the LSP case, the X-ray emission is flat and due to the 
rising part of the Inverse Compton component. At 
these relatively low energies it occurs in the Thom- 
son regime; in the ISP case, the X-ray band includes 
both the tail of the synchrotron emission and the rise 
of the Inverse Compton component. Finally, in the 
HSP case, where the emitting particles are acceler- 
ated at much higher energies than in LSPs, the syn- 
chrotron emission dominates the observed flux in the 
X-ray band and the inverse Compton scattering oc- 
curs in the Klein-Nishina regime (see FiglT]). 
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FIG. 1: Definition of different blazar types based on the 
peak of the synchrotron component (vpeak) in their SED. 
Low Synchrotron Peaked blazars, or LSP are those where 
i'^peak) is located at frequencies lower than 10^* Hz (lower 
dotted line) , Intermediated Synchrotron Peaked sources or 
ISP 10" < Upeak <10^^ Hz (SEDs with peak within the 
gray area). 



III. RADIO TO 7-RAY SEDS 

The 48 sources chosen from the LBAS list to study 
their broad band SED using simultaneous or quasi 
simultaneous radio to high-energy 7-ray data, were 
chosed on the basis of the availability of Swift observa- 
tions carried out between May 2008 and January 2009 
(which have been scheduled largely independently of 
Fermi results) and not on brightness level or any other 
condition that could influence the shape of the SED. 
This was checked by verifying that the distributions 
of redshift, optical. X-ray and 7-ray fluxes are all con- 
sistent with being the same in both subsamples. The 
quasi-simultaneous SEDs have been built in the usual 
Logi^-Log(t/F(i^)) representation. A detailed descrip- 
tion of the data analysis performed to derive the SEDs 
is reported in It is worth to underline that all 
SEDs exhibit the typical double hump distribution, 
with the first bump attributed to synchrotron radi- 
ation and the second one likely due to one or more 
components related to inverse Compton emission. 

FigUl and [3] show the SEDs of two sources among 
the 48 LBAS: 3C66A and MRK501 which can be clas- 
sified as an ISP and HSP respectively. A comparison 
between the two SEDs shows the differences in the 
position of the peaks of the two components and of 
the relative peak intensities. This behaviour is gen- 
erally found in all SEDs. Large variability is also 
present, especially at optical/UV and X-ray frequen- 
cies. Gamma-ray variability cannot be evaluated as 
the Fermi data used for our analysis are averaged 
over the entire period of data taking. The 7-ray vari- 
ability of the LBAS blazars is separately discussed in 
[Tj, while a complete description of the 7-ray spec- 
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FIG. 2: The SED of OFGL J0222.6+4302 = 3C66A . Fermi 
data appear as black stars; the quasi-simuhaneous data at 
other wavelenghts are represented by filled black squares; 
the non-simultaneous archival measurements are shown as 
gray small dots; finally filled black triangles show AG- 
ILE data. The dashed lines represent the best fits to 
the Synchrotron and inverse Compton part of the quasi- 
simultaneous SEDS (see 0] for details). 
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FIG. 4: The aox-ciro plot of the LBAS blazars sample 
(large symbols) compared to the blazars from BZCAT cat- 
alog which have radio optical and X-ray information (small 
red symbols). 
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FIG. 3: The SED of OFGL J1653.94-3946 = MRK501 . 
Fermi data appear as black stars; the quasi-simultaneous 
data at other wavelenghts are represented by filled black 
squares; the non-simultaneous archival measurements are 
shown as gray small dots; finally filled black triangles show 
AGILE data. The dashed lines represent the best fits to 
the Synchrotron and Inverse Compton part of the quasi- 
simultaneous SEDS (see for details). 



tral shape of the LBAS sources is given in [H . From 
the results of the latter paper, we note that often 
Fermi data cannot be fit by a simple power law as 
significant curvature is detected. Convex curvature is 
often observed in sources where synchrotron peak is 
located at low energies (e.g. PKS0454-234, PKS1454- 
354 and PKS1502-H106, 3C454.3 etc) whereas very flat 
or even concave type curvature is exhibited by high 
synchrotron peaked objects (e.g. 3C66A, see Fig[21 
PKS0447-439, lES 0502H-675 and PG 1246-^586 ). 



IV. OBSERVATIONAL PARAMETERS OF 
BLAZAR SEDS 

The key observational parameters that character- 
ize the SEDs of the LBAS sources have been esti- 
mated: the radio spectral index (a,), the peak fre- 
quency (t'peaj,) and the peak intensity flux of the 
synchrotron component (t'peafc-'^C^peafc))' ^'^^ peak fre- 
quency (t'p^fc) and the peak intensity flux of the In- 
verse Compton component (t'p£feF(i/p£j,)) 

The radio to microwave spectral slope has been 
measured by performing a linear regression of all ra- 
dio flux measurements used for the SEDs, both quasi- 
simultaneous and historical ones. The results of this 
evaluation is that the measured slope in our SEDs is 
quite flat (< >'^0) and consistent with being the 
same in all blazar types. 

FigH] shows the plot of the optical to X-ray spec- 
tral slope OLox versus the radio to optical spectral 
slope oiox- The plot includes all blazars in the BZCat 
catalog for which we have radio, optical and X-ray 
measurements (small red dots). Fermi FSRQs (filled 
circles), like all FSRQs discovered in any other en- 
ergy band, are exclusively located along the top-left 
/ bottom-right band, whereas BL Lacs (open circles) 
can be found in all parts of the plane, even if with a 
prevalence in the horizontal area defined by values of 
aro between 0.2 and 0.4, which is where HSP sources 
are located Q. All 7-ray selected blazars are located 
in regions covered by previously known blazars. No 
new 7-ray type of blazars has been found, in particular 
there is no evidence for the hypothetical population of 
Ultra-High energy peaked blazars (UHBLs), with syn- 
chrotron peak in the 7-ray band {i^p^^i^ >10^"). 
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A. Synchrontron peak determination 



The peak frequency {i^p^ak) ^^'^ peak intensity 



peak 



F(^'peafe) of synchrotron component of the 



SEDs have been estimated by fitting the synchrotron 
dominated part of the SED (from radio to UV in LPS 
and from radio to X-ray band in HSP blazars) with a 
simple third degree polynomial function: 



i^F^ = a ■ + h ■ u + c - V + d 



(1) 



The synchrotron peak frequency {v^ 



peak 



and the rel- 



ative peak intensity v^^ak^ i^peak) '^^'^ be alternatively 
evaluated through an empirical method that makes 
use of the spectral slopes Oiro between 5 GHz and 
5000 A and aox between 5000 A and 1 kcV with aab 
defined as: 



Q-ab 



(2) 



where fa is the rest-frame flux at frequency Va ■ The 
value of fpgafe i'^ ^^le SED of a blazar can be analyti- 
cally determined from the following relationship: 



Logivpeak) 



13.85 + 2. 30X if X<0 and Y<0.3; 
13.15 + 6. 58y otherwise. 

(3) 

-0.155-a„^ and Y^l.O- 



where X=0.565-1.433-a^o- 

0.661-Q;ro-0.339-Q;o3;. 

This relationship has been calibrated using the 



peak 



values measured from the 48 quasi-simultaneous 
SEDs and the corresponding aox and a^o values. 

The results are reported in FiglS] Despite the fact 
that aox and aox are based on non-simultaneous liter- 
ature data, the distribution of the diff'erence between 
the values obtained with the two methods has a mean 
value of 0.04 and a standard deviation of 0.58, imply- 
ing that the value of Log(i'pg^j,) can be derived even 
from non-simultaneous values of aox and aro within 
0.6 decade at one sigma level and within 1 decade in 
almost all cases. It is worth to underline that this 
method assumes that the optical and X-ray fluxes are 
not contaminated by thermal emission from the disk 
or accretion. In blazars where thermal flux compo- 
nents are not negligible, the method described above 
may lead to a significant overestimation of the posi- 
tion of i/^^^fc. 

The peak flux i^peak^i'^peak) '^^^ b^ estimated using 
the following relationship 



Log{v^F{v^) = Q.b-Log{v^)-2QA+Q.<d-Log{RbGHz) 

(4) 

where R5GHz is the radio flux density at 5GIIz in 
units of mJy. FigE] reports the value of 




Log (BsEimalBd from SEDs] 



FIG. 5: The synchrotron peak frequency Vpeak value esti- 
mated using the method based on cxox and cxro values is 
plotted against the value estimated from the SEDs. 
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FIG. 6: The synchrotron peak intensity 



!lfcF(i/. 



peak V peak) 



peak } 

value estimated using the method based on otox and Uro 
values is plotted against the value estimated from the 
SEDs. 



estimated with the two methods. Also in this case, 
the scatter around the solid line, representing perfect 
match, has an average value of -0.01 for the difference 
between the two estimates and a standard deviation 
of 0.33. It is quite remarkable that one can derive the 
synchrotron peak flux simply from Vp^^^i^ and from the 
radio flux as this implies that within a factor of 10 the 
radio emission represents a long-term calorimeter for 
the whole jet activity and the basic source power. 



B. Inverse Compton peak determination 



The peak frequency v^'^^^. and the corresponding 

peak intensity I'peak^ i^peak) inverse Compton 

bump have been also estimated by fltting the X-ray to 
7-ray part of the SED dominated by inverse Compton 
emission using the polynomial function of equation 
IIV Al However, for some sources the soft X-ray band 
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FIG. 7: The 7-ray power law photon index (F) plotted 
against the log of synchrotron peak energy. A clear corre- 
lation is present. 
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FIG. 8: The 7-ray power law photon index (F) plotted 
against the log of Inverse Compton peak energy. A clear 
correlation is present. 



is still dominated by synchrotron radiation, and only 
the Fermi data can be used to constrain the inverse 
Compton component and the above method is sub- 
ject to large uncertainties. For this reason, in these 
cases, the simultaneous data points have been fitted to 
a SSC model with a log-parabolic electron spectrum 
[| using the ASDC SED [H interface to fit. 

The i^peak^i'^peak) bccn determined as the value 
of V that maximizes the vF^v) in the IC polynomial 
function or the predictions of the SSC model. The 
best fit to both the synchrotron and inverse Compton 
components appear as dashed lines in Fig. [2]and[3l 

FigE and [8] respectively show that both the i^p^^k 
and the t'p^j. present a strong correlation with their 
7-ray spectral slope F taken from From the two 
plots is clear that the scatter is larger in the regions 

of low Ji'c'' -steep values of F, probably reflecting the 
presence of 7-ray spectral slope curvature. The best 
fit to the I'p^ak'^ relationship is 

Log{,y'^Z,) = -4.0-r + 31.Q (5) 

Since the 48 sources for which we have quasi- 
simultaneous SEDs are representative of the entire 
LBAS sample, the above equation has been also used 
to estimate the of the LBAS sources for which we 
have no simultaneous SEDs [5i] . The statistical uncer- 
tainty associated to t'p^fc calculated via ea llVBI can 
be estimated from the distribution of the difference 
between vl'-^,, measured from the SED and that from 
ea llVBI This distribution is centered on the value of 
and has a sigma of 0.51; considering that the value of 
^peak from the SED is also subject to a a similar error 
we conservatively conclude that Log(i'p^j,) estimated 
through eq llVBI has an associated error of about 0.7. 



V. DISTRIBUTION OF SYNCHROTRON 
AND INVERSE COMPTON PEAK 
FREQUENCIES 

The new classification scheme of blazars introduced 
in section |TT] and the method used to estimate fpg^fc > 
described in section IIV Al allow to classify the LBAS 
sample in terms of its content of LSP, ISP and HSP ob- 
jects and to compare it with that of samples selected 
in other energy bands. The distribution of the syn- 
chrotron peak frequency {t^pf,ak) of LBAS blazars (es- 
timated using the aox-c^ro method) shows that while 
the i/^g^j, distribution of FSRQs starts at ~10^^-^ Hz, 

peaks at lO^^ '^ Hz and it does not extend beyond 
j^-|^q14.5 jj^j the distribution of BL Lacs is much flat- 
ter, ~10, it starts at ^lO^'^ Hz and reaches much 
higher frequencies (wlO^^ Hz) than that of FSRQs. 
This result has been compared with the distribution 
of fpeafe of the sample of FSRQs and BL Lacs detected 
as foreground sources in the WMAP 3-year microwave 
anisotropy maps and with the I'pgjjj. distribution 
of the LBAS sample with that of the X-ray selected 
sample of blazars detected in the Einstein Extended 
Medium Sensitivity Survey [ll|- It can be noted that 
the t'peafc distribution of FSRQs is consistent with be- 
ing the same in the 7-ray, radio/microwave and in the 
X-ray band. Moreover, the large majority of FSRQs 
are of the LSP type while no FSRQs of the HSP type 
have been found at any frequency. On the contrary, 
the I'pga;; distribution of BL Lac objects is very differ- 
ent in the three energy bands. It is strongly peaked at 
^-[^q13.3 in the microwave band, where HSP sources 
are very rare, whereas in the X-ray and and 7-ray 
bands HSP sources are more abundant than LSPs. On 
the other side, the distribution of the inverse Comp- 
ton peak frequency I'p^g.k' of FSRQ (dot-dashed 
histogram) and the BL Lacs (solid histogram) in the 
LBAS sample. The two distributions are quite differ- 



eConf C091122 



6 



2009 Fermi Symposium, Washington, D.C., Nov. 2-5 




FIG. 9: The distribution of synchrotron peak energy for 
the sample of LBAS FSRQ (top panel) and BLLacs (bot- 
tom panel). 
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FIG. 10: The distribution of the inverse Compton peak 
energy for the sample of LBAS FSRQ (dot-dashed line) 
and BLLacs (solid line). 



ent with the BL Lacs exhibiting much higher I'p^j. 
values. This is most likely due to the same reason 
that causes the different v^„„ ,. distributions in the two 



blazar subclasses. 



peak 



VI. CONCLUSIONS 

A sample of LBAS blazars detected by the Fermi 
LAT in the first three months of scientific opera- 



tions has been selected to build simultaneous or quasi- 
simultaneous SEDs using data obtained also from 
Swift, radio/mm telescopes and infra-red and optical 
facilities. The collection of SEDs allowed us to es- 
timate important parameters characterizing the SED 
of 7-ray selected blazars, such as the frequency and 
intensity of both the synchrotron and inverse Comp- 
ton peaks, both directly from the simultaneous data 
or indirectly using a refined version of the method 
proposed in @ based on the position in the aox-oiro 
plane for the former and on the slope of the 7-ray 
spectrum for the latter, as the 7-ray spectral slope 
and v^'^^f. are strongly correlated. Among the sev- 
eral results we obtained, we can conclude that 7-ray 
selected blazars have broad-band spectral properties 
similar to those of radio and X-ray discovered blazars 
implying that they are all drawn from the same under- 
lying population; the distribution of the synchrotron 
peak frequency is very different for the FSRQ and BL 
Lac subsamples. The results rule out the existence 
of FSRQs of the HSP type, consistent with what also 
observed in radio, microwave and X-ray surveys; the 
BL Lac minimum 1/^^^^^, appears to be larger than in 
FSRQs and this could be due to some intrinsic differ- 
ence in the mechanism of particle acceleration in the 
two types of blazars or to a mere selection effect. 
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